Resonant electron transmission through a finite quantum spin chain 
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Electron transport in a finite one dimensional quantum spin chain (with ferromagnetic exchange) 
is studied within an s — d exchange Hamiltonian. Spin transfer coefficients strongly depend on the 
sign of the s — d exchange constant. For a ferromagnetic coupling, they exhibit a novel resonant 
pattern, reflecting the salient features of the combined electron-spin system. Spin-flip processes are 
inelastic and feasible at flnite voltage or at finite temperature. 



Motivation and scope: There is a growing interest in 
electronic devices which transport electron spin ^ (to- 
gether with its charge). Beside the Kondo physics 0, 
electron transmission through a region containing a lat- 
tice of magnetic atoms (mainly a domain wall) is a focus 
of experimental § g |l § and theoretical § || |] 
investigations. In this context, the basic ingredient is an 
exchange interaction gs ■ S„ between the electron spin 
operator s and an atomic spin S„ localized at a point 
Xn, with an "s-d" coupling g. Most theoretical works 
treat the atomic spins either purely classically (as local- 
ized magnetic fields) , or semi classically (the operator S„ 
is replaced by its expectation value < S„ > in the ground 
state of the quantum spin system). It is expected, how- 
ever, that the quantum nature of the magnetic system 
will be crucial at low temperatures. 
In the present work, transmission of electrons through a 
finite one dimensional quantum spin 1/2 chain (with fer- 
romagnetic exchange) is studied, starting from an s — d 
exchange Hamiltonian. The many body Kondo problem 
is avoided by an application of a magnetic field which 
removes the degeneracy of the ground state. It also en- 
ables cutting the spectrum of the spin system off at the 
one magnon level. The ensuing scattering formalism can 
then be handled within a transfer matrix algorithm, in 
terms of which the spin transfer coefficients (STC) are 
calculated as function of the electron Fermi energy. It 
is found that for ferromagnetic s — d exchange, the STC 
display a novel resonant structure which manifests the 
richness of the combined electron spin chain system. 
Hamiltonian: Consider electrons (mass m and charge 
— e) moving in a one dimensional wire (along x) inter- 
acting with a chain of quantum spins {e.g. magnetic 
atoms) S„, {n = l,2,...N, Sjv+i = Si) localized at 
points Xn = na {n = 1, 2, ...N), subject to a perpendicu- 
lar magnetic field H = HzZ. Henceforth a and ?i^/2ma^ 
are exploited such that all length and energy quanti- 
ties are dimensionless. The Hamiltonian of the system 
H = He + Hs~d + Hs contains an electronic part H,,, an 
electron-spin s — d exchange interaction Hs-a and a term 
Hs controlling the isolated spin system. The latter is, 

Hs = — ^( JS„ • S„+i + EzSnz), (1) 
n 

where J > and Ez is the atomic Zeeman energy. 



The spin Hamiltonian Hs is projected on the subspace 
spanned by its + 1 lowest energy states \k > [k = 
0, 1, . . .N). If the magnetic field is strong enough such 
that Ez > JS then, beside the (non degenerate) ground 
state |fc = > (with all spins up along z and energy 
Eq = -N{S^J + EzS) and Y.n ^nz = f ), these are the 
one magnon spin waves (OMSW) with ^„ Snz — ^ 



In terms of states \n >= S„ 
corresponding energies are, 
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Cutting the spectrum off at the OMSW states level is 
justified since above it there is a gap of 2SEz until the 
lowest energy of two magnon spin waves. 
The electronic and s — d parts of the Hamiltonian should, 
in principle, manifest the many-body aspects of the prob- 
lem (in the Kondo sense). In first quantization that 
amounts to, H^. = — — Siz and an exchange 

part Hs-d = g X^i X^rt ^i^i ~ ■ (^e is the electronic 
Zeeman energy). The sum over i runs, principally, on all 
the electrons in the wire. However, due to the presence 
of an external magnetic field, the ground state of the 
spin system is non degenerate and the Kondo effect is 
absent. Hence, the dynamics of the system is adequately 
described by the single electron s — d exchange Hamilto- 



H 



dx'^ 



heSz -\- g^^S{x ~ n)s ■ Sn + Hs 



The Schrodinger equation is then, 

H^ix;{Sn}) = E-^ix;{Sn}), 



(3) 



(4) 



where ^(a;; {S*,!}) and E are the total wave function and 
total energy (electron and spin system). 
The scattering problem: Between spins and outside 
the spin system, electrons propagate as plane waves with 
momenta po-fc = ~ Ek+ crh,,), (some of which might 

be purely imaginary), where a = ±i, (alternatively 
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cr =ti) is the electron spin projection along z. Solution of 
the scattering problem means the evaluation of transmis- 
sion amplitudes t^ik'ak{E), in which ak are initial elec- 
tron and spin system quantum numbers, while a'k' are 
the final ones. To extract the transmission amplitudes 
from the Schrodinger equation, note that the presence of 
spins at isolated points suggests using a transfer matrix 
formalism. For Xn-i < x < Xn the total wave function is 
expanded in electron spinors Xa and spin states |fc >, 



(5) 



Analogous expansion holds for Xn < x < Xn+i- The 
matching conditions at x = n read, 



- n~;{Sn}) = ^{x = n+; {5„}), 



(7) 



^'{x - n+; {Sn}) - ^'{x = n-;{Sn}) 

= gs-S„*(a; = n-;{5„}). (8) 

Employing the expansion (|^) and applying on the left 
the bra < fc| (g) xJ-j the relations (0J|) imply 4(iV + 1) 
equations from which the coefficients an+iak,bn+iak arc 
expressible in terms of an^k , bnak ■ At the end of this pro- 
cedure (which is briefly detailed here) , it yields the trans- 
fer matrix t„ that carries the system from n — to n -I- 0. 
To be definite, the order of coefficients (which determines 
the structure of transfer matrices) is (a„f , a„| , b„| , b„|). 
The required manipulation is to evaluate matrix elements 
of four (iV-l-l)x (A^-|-l) matrices la-'a defined as operators 
in spin wave space, 



< /c'|/av|fc >=< k'\ (g) xl'S ■ S„Xa (g |fc > 



(9) 



Straightforward calculations yield, (separating fc, /c' = 
and fc, fc' ^ 0) 



< fc'l/nl^ 



<k'\Iu\k>-- 
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The transfer matrix across x = n (t„) can now be writ- 
ten down, remembering that in many channel problems, 
the plane waves are normalized to have unit velocity. In 
terms of the 4(A^ + 1) x 4(A^ + 1) diagonal matrix. 



p 2 = diag{p^ ',P| ',Pt ',Px 



(11) 



where Po- = {Pa-OjPai, --PaN) is the vector of -I- 1 chan- 
nel momenta, the result is. 
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Propagation between spins is controlled by the diagonal 
4:{N -I- 1) X 4(A^ -t- 1) matrix of phases, 



ipi 



)■ 



(13) 



The transfer matrix T„ across a unit cell (from ti — 1 + to 
n + 0) and the total transfer matrix T are then, r„ — Ar„ 
and T = 11^^=1 ^n- With some modifications required 
in case there are evanescent modes, it is not difficult 
to show that in the above 4 block partition of trans- 
fer matrices into a and b sectors, the transfer matrices 
T„, Tn and T satisfy the current conservation constraint, 
TnO'zT^i = Tz- Finally, the matrix of amplitudes tcr'a with 
elements ta-'k'ak{E) is simply given by T~J^. 
Conductance and spin transfer coefRcients: An ex- 
perimentally relevant question is the following: Applying 
a small voltage V across the wire between and A^ + 1 
and letting a unit flux of electrons at Fermi energy ep 
and spin component a reach it from the left, what is 
the current la-' a of electrons with spin a' in the sys- 
tem? The ratio Ga'a-isp) = la'a/V is the corresponding 
STC. The conductance per spin is G(j — ^^a' Ga'a and 
the conductance is G = Scr^o'- zero temperature 
(and for infinitesimal voltage drop V) the system is in its 
ground state and only elastic scattering (of the electrons) 
is allowed. Due to spin conservation, processes involving 
spin fiip are inelastic and hence are forbidden. Thus, 
the quantities to be calculated at zero temperature and 
infinitesimal voltage difference are, 



(14) 



where ep + Eq = E is the total energy of the sys- 
tem. Evidently, G|| is trivial since it involves only 
the ground state |0 > of the spin system. Coupling 
terms to OMSW, < fc ^ 0| ® xl'S ■ S„XT «> |0 >= 
since the initial state X| (8) |0 > is stretched with to- 
tal spin N/2 + 1/2 while the final state is not. There- 
fore, the quantum spins are replaced by their averages 

< 0|S„|0 >= Snz and spin up incoming electrons are 
scattered by a sequence of A^ identical delta function po- 
tentials of strength 17/4. On the other hand, in- 
volves virtual OMSW excitations and hence it is non 
trivial. For examples, it gets second order contributions 

< 0| (g) xjs • SmXi ® 1^ >< ^1 ® xjs • S„xx (g) |0 >. 

Of special interest are spin fiip processes. If the po- 
tential drop eV exceeds the gap between the ground 
state energy Eq and the lowest OMSW energy (that is. 



2 



eV > 2SEz), then the process T<~i is energetically fea- 
sible even at T = 0. The electron stays at the Fermi 
energy and the spin system is heated from Eq to Ei , the 
energy being supplied from the external voltage. Assum- 
ing that linear response is still valid, the corresponding 
STC is G|x(eF) |ttixo(£F + ^'o)^- After the system 
is excited from the ground state to the lowest OMSW 
state, the process t^i is blocked. It can occur again 
only if the spin system returns to its ground state. This 
can happen for example by the process i<— t (which, by 
spin conservation, is allowed if the spin system is in a 
OMSW state but is otherwise forbidden on the ground 
state). Finally, consider the spin- flip process encoded in 
^ni^p) = \tioTk{£F + Ek)\'^. This can happen only at 
T > where the spin system contains also a OMSW com- 
ponent. At the end, the spin system cools down to the 
ground state, but regains the excited component energy 
from the heat bath. 

Possible experimental realization: Out of the four 
parameters g, J, H and ep, the last two are the ones 
which are easiest to handle and control. However, the 
strength of the magnetic field is somewhat constrained 
by the condition Ez > JS justifying projection on the 
OMSW subspace. The s — d coupling g is material depen- 
dent and can have either sign but in general \g\ ^ | J|. 
The atomic and electronic Zeeman energies are similar in 
magnitude, Ez ^ h^. Thus, one should concieve a ma- 
terial consisting of itinerant electrons and localized mag- 
netic moments in which the Fermi energy can be tuned 
by a sutable gate. Possible candidates are magnetic semi- 
conductors {e.g. FeSi) or strongly correlated semimetals 
such as CeP and CeSb. 

Results: Calculations are performed on a chain of 
= 12 spins, using JS — 0.06, and taking g — ±1 
for the electron-atom s — d exchange strength. The 
atomic and electronic Zeeman energies are taken to be 
Ez = he = 0.2. The Fermi energy is taken with refer- 
ence to the ground state of the spin system in units of 
\J\S. The results at T = are presented in figures 0, ||, 
and H below. 
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FIG. 1. Non-flip spin transfer coefflcient G]] {eF) (equa- 
tion p^ ) for electron transmission through a chain consisting 
of A*' = 12 spins prepared at the ground state |0 >. Solid 
(dashed) lines correspond to g = —1, a, ferromagnetic (<? = 1, 
an antiferromagnetic) s — d exchange coupling. Other param- 
eters are J = 0.12, Ez = he = 0.2. 
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FIG. 2. Non-flip spin transfer coeflicient Gxx(eF) (solid 
line) and spin flip transfer coeflicient G|x(£f) (dashed Hne) 
for a ferromagnetic s — d exchange coupling g = —1. Other 
specifications are as in figure hi 



0.14 




FIG. 3. Same as figure g for g = 1 (an anti ferromagnetic 
s — d exchange coupling). 

As asserted above, the behavior of G^^ (figure |l|) is 
somewhat expected since only the SzSnz interaction is 
effective. For a ferromagnetic s — d coupling the electron 
encounters a periodic potential composed of attractive 
S functions while for an antiferromagnetic exchange the 
potentials are repulsive. This leads to a band structure 
determined by, \cosko + -^sinkgl < 1, where ko = y/sp- 
It is easy to see that the above range oi ep is inside the 
band. 

The patterns of the STC Gn and G-fx (figures ^ for 
g = —1 and |^ for g = 1) have a much richer content. 
In particular, for a ferromagnetic exchange coupling (fig- 
ure ^ they display a series of very narrow resonances. It 
is quite remarkable that in a magnetic field and with fer- 
romagnetic s — d coupling, still, G^ approaches unity at 
resonance. These resonances are novel, that is, they are 
not the usual ones encountered in resonance tunneling, 
since there is neither potential scattering here nor dou- 
ble barrier. As discussed above, the pertinent processes 
involve virtual excitations of the OMSW, and hence, it 
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reflects the structure of the electron spin chain system. 
Unlike potential scattering in which the scattering is ef- 
fected upon static (quenched) impurities, we encounter 
here scattering from a quantum system with an inter- 
nal structure and many non-degenerate levels (the case 
of degenerate levels is exemplified by the Kondo effect). 
It is also noticed that these resonances appear only for 
a ferromagnetic s — d coupling g < 0. Some light can 
be shed on these points by inspecting the corresponding 
bound state problem for a closed system. Namely, the 
one dimensional wire is closed into a ring of length N 
and the bound state energies (pertaining to the total sys- 
tem, electron and spin chain) can be computed by solving 
the equation det[T{E) — 1] = (in order to avoid orbital 
effects, the magnetic field is tuned to yield an nteger num- 
ber of flux quanta through the ring). While there is no 
simple relation between an open wire conductance rson- 
ances and bound states on a closed wire, it is nevertheless 
instructive to notice that the bound state energies in the 
ferromagnetic case are lower than those in the antiferro- 
magnetic case. In fact, the latters fall mostly outside the 
range of Fermi energies examined in these simulations. 
At finite temperature, the spin system can be activated 
to the lowest OMSW (which, unlike higher OMSW, is 
not degenerate). The process | 1 is then feasible 
and the relevant STC Gj-f is depicted in figure 
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FIG. 4. The spin transfer coefficient (ef) for a spin sys- 
tem prepared at the first OMSW state. The value of the s — d 
exchange constant is g = =pl (solid and dashed lines). Other 
parameters are as in figure ^. 

It mainly reflects the coupling between the ground and 
lowest OMSW state of the spin system. Higher OMSW 
are virtually excited in higher orders through non-fiip 
processes, e.g. (i 0) ^ (t /c > 1) <— (t 1). The resonant 
pattern prevails also in this case. Note that the actual 
value of the STC should be reduced by the activation 
factor e-'3(-Ei--Eo). 

Conclusions: The formalism for studying electron 
transport in a quantum spin chain is developed. Appli- 
cation of an external magnetic field removes the Kondo 
problem as it hfts the degeneracy of spin states. (For 
infinite chains it also restores long range order which is 



otherwise absent due to the Mermin - Wagner theorem). 
Spin transfer coefficients G^'a are evaluated and analyzed 
as a function of the Fermi energy. Non-flip processes 
t^t and i^i are elastic. The structure of is trivial 
since the problem is equivalent to the scattering of polar- 
ized electrons on a sequence of delta function potentials. 
On the other hand, the behavior of is much richer. 
For a ferromagnetic s — d exchange coupling it displays 
a novel resonant pattern which is distinct from the one 
encountered in potential or double barrier transmission. 
That is, it reflects the complex structure of the electron 
spin-chain system. Spin-fiip processes are inelastic. The 
transition t^i is endothermic, and is feasible at T = if 
the potential difference across the spin system overcomes 
the energy gap between the ground and lowest OMSW 
states. The transition i<— t is exothermic and feasible at 
T > if OMSW are activated. Although spin flip coef- 
ficients are small relative to non spin fiip ones, they are 
still sizable. In passing, we note that there are interest- 
ing theoretical works discussing the pertinent many body 
problem @ 
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